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Abstract 

We study the numerical differentiation formulae for functions given in grids with arbitrary 
number of nodes. We investigate the case of the infinite number of points in the formulae for 
the calculation of the first and the second derivatives. The spectra of the corresponding weight 
coefficients sequences are obtained. We examine the first derivative calculation of a function 
given in odd-number points and analyze the spectra of the weight coefficients sequences in 
the cases of both finite and infinite number of nodes. We derive the one-sided approximation 
for the first derivative and examine its spectral properties. 
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Introduction 

The problem of numerical differentiation is a long-standing issue. There are plenty of published 
works devoted to the generation of finite difference formulae in both one and multi dimensional 
lattices (see, e.g., Ref. [I]). However, many of those methods require preliminary construction of 
an interpolating polynomial, and hence are very awkward. Moreover, the majority of the previous 
techniques are valid in the case of a function given in the limited number of nodes. 

The finite difference formulae for the calculation of any order derivative in a one dimensional 
grid with arbitrary spacing were discussed in Refs. However, only recursion relations for the 

weight coefficients have been derived. The explicit formulas for the derivatives calculation were 
recently derived in Ref. [4] on the basis of the generalized Vandermonde determinant. 

It should be noted that the low order derivatives (the first and the second ones) as well as 
equidistant lattices are of the major importance in many problems of applied mathematics and 
physics. The first and the second numerical derivatives in the equidistant one dimensional grid 
were studied in Ref. [5 J . The finite difference formulae for the central derivatives of a function given 
on a lattice with arbitrary number of elements have been derived in that work. It is important 
that these formulae have been obtained in the explicit form. This method enabled one to examine 
the spectral properties of weight coefficients sequences as well as to analyze the accuracy of the 
numerical differentiation. 

In the present paper we continue to study the numerical differentiation formulae for functions 
given in grids with arbitrary number of nodes. On the basis of the results of Ref. [5] in Sec. [1] 
we investigate the case of the infinite number of points in the formulae for the calculation of the 
first and the second derivatives. The spectra of the corresponding weight coefficients sequences 
are also obtained. Then, in Sec. [2] we examine the first derivative calculation of a function given in 
odd-number points. We also analyze the spectra of the weight coefficients sequences in the cases 
of both finite and infinite number of nodes. In Sec. [3] we derive the one-sided approximation for 
the first derivative and examine its spectral properties. It is worth noticing that the derivations of 
the finite difference formulae in all cases are performed for the arbitrary number of points. Finally 
in Sec. 2] we resume our results. 
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1 Spectral properties of the first and the second derivatives 
for infinite number of points 



Let us study the function f(x) given in the equidistant points x m , f(x m ) — f m , where to 
I ) ±n. It was found in Ref. [5j that the first and the second derivatives are approximated as 
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where /i is the distance between nodes. The coefficients aS'(n) and ctm'(u) can be calculated 
explicitly for arbitrary n (see Ref. [5]). 

The spectral properties of the sequences o<m\n) and a„ (n) in Eqs. fll.ljl and (|1.2[) in the case 
of finite number of interpolation points were carefully examined in Ref. [5]. We found out that 
the more points we involved in the sequence (n) the more close to linear the corresponding 
spectrum was. Thus the considered sequence produces more accurate first derivative of a function 
in the case of great number of points. As for the sequence a m (n), it was also shown in Ref. [5] 
that its spectrum approached to parabola if n > 1. We expect that the corresponding spectra will 
be exactly linear and parabolic ones if n — > oo. 

Let us consider the spectral characteristics of the sequences a^i (n) and a>m {n) in the case 
of infinite number of interpolation points. First we remind the result for the a^(n) in the limit 
n — > oo (see Ref. [5]) 

a®= lim ag\n) = (-l) m+l -. (1.3) 
The Fourier transform of a function f(x) can be presented in the form (see, e.g., Ref. [5]) 

c(w) =hJ2 e ~ iux f( x ) = h E e- iumh f(mh). (1.4) 
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The inverse Fourier transformation is given by the following expression: 
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and has the cutoff at high frequencies, |w| < ir/h 
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Now we can calculate the spectrum of the sequence am , 
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where we use Eqs. (|1.3j) and (|1.4|) . Note that Eq. (|1.5[) is valid if < a; < 7r//i. The first derivative 
of the function can be expressed via the spectra /3i(w) and c(u>), 

f(x) = k fS h ^h pi{uo)c{ujyuix > x = kH - (i - 6) 

Using the result for the calculation of f3\(uj) presented in Eq. (jl.5| we readily find that 

/'(*)= r 1 ' 1 ^(iu)c{u)e™*. (1.7) 
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Figure 1: The spectra of differentiating niters, (a) a m and (b) a m , in the case of infinite number 
of points. 



Eq. (|1.7p shows that the first derivative calculation with help of the sequence a m gives the exact 
value of the derivative in the case of infinite number of interpolation points for all frequencies 
except w max = w/h. The fact that the first derivative computation does not give correct results at 
u> = u> max also follows from Fig. [IJa). However, it can be verified directly with help of Eq. (|1.3p 
for the function f m = (—1)"' = cos(w max m/i). 

In Ref. [5] we showed that the use of the coefficients a™ give exact value for the first derivative 
of the function y(x) = sin(w ma x^/2). However, Eq. (|1.7|) [see also Fig. [T^a)] indicates that this 
method will give the correct results not only for uj = u; max /2 but also for all frequencies to < w max . 

The second derivative calculation in the case of infinite number of interpolation points can be 
analyzed in the similar manner as we have done it for the first derivative. The explicit form of the 
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sequence ct m is 
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For the spectrum of the sequence a m we obtain 
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It should be noted that Eq. (|f .8|) is valid for all frequencies < w < ir/h. The expression for the 
second derivative takes the form 

n*) = i r /h i^h - %m ^y ux = C' h ^(-^ 2 ) C (^)e-, x = kh. (i.9) 

Eq. (|1.9p demonstrates that the computation of the second derivative with the use of the sequence 

(2) 

am gives the exact results in the case of infinite number of interpolation points for all frequencies 
even including the maximal one. The spectrum 02(u>) is depicted in Fig. [TJb). 

2 The first derivative computation of a function given in 
odd-number points 

In this section we discuss the calculation of the first derivative in the case of a function given 
in odd- number nodes. Then we discuss the spectral properties of the derived weight coefficients 
sequences in the case of both finite and infinite number of nodes. 
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Figure 2: (a) The spectra of various sequences aim (n) at N — 2000. (b) The example of the 
use of the sequence a£ (n). A rapidly oscillating function corresponds to the dashed lines, its 
envelope functions are shown by the dotted lines. 



It follows from Fig.flja) that the computation of the first derivative gives unsatisfactory results 
at high frequencies near w max - In order to introduce the numerical differentiation of such rapidly 
oscillating functions we consider the modified sequence 
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4m+iO) = — 77 (i/2), I ■ m = 0,...,n-l, (2.1) 
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and Oi^^in) = 0. 

The coefficients in Eq. (|2.ip can be formally derived if we consider the first derivative calculation 
of a function given in the odd-number points only 
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Note that originally the function f(x) was given in 2n + 1 points. 

It is worth noticing that the results for the computation of the weights with help of Eq. (|2.1| 
in some particular cases (namely for n = 3, 5, 7 and 9) coincide with those presented in Ref. [2] for 
the centered approximations at a 'half-way' point. However, the method for the central derivatives 
calculation elaborated in our paper enables one to get the expressions for the weight coefficients 
in the explicit form for any number of nodes. 

We consider the spectral properties of the obtained sequence a>m (n). Using the technique 
developed in Ref. [5] one can compute the spectrum of the sequence in question, 
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The spectra of the sequences o?m^ (n) are depicted in Fig. [21(a) for the various values of n at 
N = 2000. The function yi/ 2 (r) has the form 
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Figure 3: The spectrum of the sequence am in the case of infinite number of points. 

_ j2irr/N, if < r < AT/4, 

2/1/2 W - |^ 27rr/iV, if iV/4 < r < AT/2. 

It follows from this figure that for n — 1 the imaginary part of the spectrum is the function 
sin(27rr/iV). The linearity condition is satisfied only in the vicinity of zero and N/2. However at 
n = 10 linearity condition remains valid for r < 350 and r > 650. 

We examine the details of the differentiation procedure performed with help of the sequence 
ani 2 \n). If a function is slowly varying, then Eq. f|2 . 2|) gives the approximate value of the first 
derivative. It also results from Fig. [UJa). However, if a function is rapidly oscillating (e.g., at 
lu < w max ), we can consider its upper and lower envelope functions. Therefore the use of the 

sequence a™^ 2 ' (n) will produce the first derivative of envelope functions since in Eq. (|2.2|) we 
calculate the sum in odd- number points only. Envelope functions can be treated as "smooth" in 
this case. Fig. HJb) schematically illustrates this process. 

Now let us discuss the case of infinite number of the interpolation points. We can treat the 
sequence al^/ 2 ^ (n) in the similar way as it was done in Ref . [5^ . Indeed, proceeding to the limit 
n — > oo in Eq. (|2.1[) we find that 
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In Eq. (|2.3|) we used the known value of infinite product, 
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With help of Eq. (J2T3J) it is possible to obtain the spectrum of the sequence a 
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if < lu < n/2h, 
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The imaginary part of the spectrum /3i/a(w) is presented in Fig. [3] As it follows from Eq. (|2.4p (see 
also Fig. [3]) the sequence cSm 2 ^ performs the differentiation of a function in question if lu < w max /2, 



or its envelope functions if lu > u> n 
Ref. 0. 



s /2. The of case lu 



</2 was also considered in details in 



3 One-sided approximation for the first derivative 



In this section we derive the weight coefficients for the one-sided approximation of the first deriva- 
tive and then we analyze the spectral characteristics of the weight coefficients sequence. It should 



be noted that the derivation of the weight coefficients is analogous to case the of the central 
derivatives which was carefully examined in Ref. 5J. 

Without restriction of generality we suppose that we approximate the first derivative in the 
zero point. Let us consider the function /(x) given in the equidistant nodes x m = mh > 0, where 
to = 0, . . . ,n, and h is the constant value. We can pass the interpolating polynomial of the nth 
power through these points, 

n 

P n{x) = ^C k X k . 
k=0 

The values of the function in the nodes x m = mh, f m = / (x TO ), should coincide with the values 
of the interpolating polynomial in these points, 
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In order to find the coefficients Ck, k = 0, . . . ,n, we receive the system of inhomogeneous linear 
equations with the given free terms f m . It will be shown below that this system has the single 
solution. 

We will seek the solution of the system (13. ip in the following way: 
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where »m (n) are the undetermined coefficients satisfying the condition, 
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It is worth to be noted that, if we set k = and '^0 in Eq. 
should be imposed on the coefficients a m ' (n) 



, we obtain the constraint which 



£>«(n) = 0, 1 = 1, 
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Analogous relation between the weight coefficients was established in Ref. [2]. In deriving of 
Eq. (|3.3p (as well as in all subsequent similar formulae) we suppose that to = 1 if to = 0. 
Let us resolve the system of equations (|3.2[) according to the Cramer's rule 
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(3.4) 



(3.5) 
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In Eq. I|3.5p we use the formula for the calculation of the Vandermonde determinant. From Eq. 
it follows that the determinant of the system of equations (|3.2p is not equal to zero, i.e. the system 
of equations (|3.ip has the single solution. 

The most simple expression for Am (n) is obtained in the case of I = 1 that corresponds to the 
calculation of the first-order derivative 
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(j-i), m=l,...,n. (3.7) 
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From Eq. (|3.4p as well as taking into account Eqs. (|3.5|) and (|3.7p we get the expression for the 
coefficients a m (n) 

M)f n \ = I i\m+l 1 ( n 
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where 
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are the binomial coefficients. It is remarkable to note that the coefficient (n) = n. To simplify 



numerical calculations (especially the analysis of the spectra of the derived sequences) Eq. (|3.8p 
should be rewritten in the form 

"m(n) = Vt' m=l,...,n, 

mp m (n) 

where 
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In order to find the coefficient <Zg (n) we use Eq. (13. 3p rather than compute the determinant 
(|3.6p . Thus we receive the following expression for this coefficient, 
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Eqs. p.8p and (|3.9p provide the weight coefficients for the one-sided approximation of the first 
derivative of the function f(x) given in n + 1 equidistant nodes, 
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The results for the computation of the weights in some particular cases (namely for n = 1, 2, . . . , 8) 
coincide with those presented in Ref. [2]. However, the technique for derivatives calculation de- 
veloped in the present work allows one to obtain the expressions for the weight coefficients in the 
explicit form for any n. 

Without derivation we mention that on the basis of Eqs. (|3.4p ~ (|3.6p one can find the coefficients 
at' (n) that correspond to the computation of the nth-order derivative 

a^(n) = (-ir+ n -( n ). m = 0,...,n. (3.10) 
n! \m J 

It should be noted that Eq. (|3.10p is consistent with Eq. 




Now we consider the spectral properties of the derived sequence dm (n). Using the results 
of the previous section (see also Ref. [5]) we readily find the expression for the spectrum of the 
considered sequence, 

h(r) = 8 m(") ex P (~ i ~^~ rnr ) ■ 

The imaginary parts, 3fm[6*(r)l , of the spectra of the sequences am\n) for the various values of 
n at N = 2000 as well as the linearly growing sequence I(r) — 2irr/N are presented in Fig. 01a). 
It follows form the this figure that the imaginary parts are close to the linear sequence only in the 
vicinity of zero (r < 200) even at n = 5. It points out that the one-sided approximation of the 
first derivative has worse accuracy in comparison with the central derivatives. This fact was also 
mentioned in Ref. [7J. Therefore, the application of the sequence ujj'(n) for the calculation of the 
one-sided first derivative will give reliable results only for slowly varying functions. 

The spectrum b\(r) has not only imaginary part, but also nonzero real part since <ig (n)^0. 
The real parts, 5Re [b* (r)] , of the spectra of the sequences a£)(n) for the various values of n at 
N = 2000 as well as the constant sequence R(r) = are shown in Fig. [3Jb). It can be seen from 
this figure that the real parts of the spectra are close to zero if r < 100 at n = 1, and if r < 300 at 
n = 3 and n — 5. The deviation from zero is especially great if r > 700 at n = 3, and if r > 500 
at n = 5. Such a behavior of the real parts of the spectra also reveals the limited level of accuracy 
of the one-sided first derivative. 



4 Conclusion 

In conclusion we note that in our paper we have studied the numerical differentiation formulae 
for functions given on grids with arbitrary number of nodes. In Sec. Q] we have investigated the 
case of the infinite number of points in the formulae for the calculation of the first and the second 
derivatives. The spectra of the corresponding weight coefficients sequences have been obtained. 
It has been revealed that the calculation of the first derivative with help of the derived formulae 
gave reliable results for all spacial frequencies except w max - As for the calculation of the second 
derivative we have shown that the corresponding formulae were valid for all spacial frequencies 
including w m ax- In Sec. [2] we have examined the first derivative calculation of a function given in 
odd-number points. We have also analyzed the spectra of the weight coefficients sequences in the 
cases of both finite and infinite number of nodes. It has been found out that the obtained formulae 



perform the differentiation of the considered function if u> < w max /2, and its envelope functions if 
u > w max /2. In Sec. [3] we have derived the one-sided approximation for the first derivative and 
examined its spectral properties. The accuracy of the one-sided first derivative has been discussed. 
On the basis of the spectral properties of the weight coefficients sequences it has been shown that 
the accuracy of the one-sided approximation for the first derivative was essentially lower compared 
to the computation of the central derivatives. This our result is in agreement with previous works 
(see, e.g. Ref. [7]). Nevertheless, the obtained one-sided first derivative formulae could be of 
use in solving differential equations by means of numerical methods. It is also possible to apply 
the elaborated technique of the numerical differentiation in construction of quantum field theory 
models of unified interactions. 
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